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Although numerical methods are required to evaluate the stransverse mass, Mti, for general input 
momenta, non-numerical methods have been proposed for some special classes of input momenta. 
One special case, considered in this note, is the so-called "massless limit" in which all four daughter 
objects (comprising one invisible particle and one visible system from each "side" of the event) have 
zero mass. This note establishes that it is possible to construct a stable and accurate implementa- 
tion for evaluating Mti based on an analytic expression valid in that massless limit. Although this 
implementation is found to have no significant speed improvements over existing evaluation strate- 
gies, it leads to an unexpected by-product: namely a secondary variable, "My™™", that is found 
to be very similar to Mt2 for much of its input-space and yet is much faster to calculate. This is 
potentially of interest for hardware applications that require fast estimation of a mass scale (or QCD 
background discriminant) based on a hypothesis of pair production - as might be required by a high 
luminosity trigger for a search for pair production of new massive states undergoing few subsequent 
decays (eg di-squark or di-slepton production). This is an application to which the contransverse 
mass Mct has previously been well suited due to its simplicity and ease of evaluation. Though 
Afy2 Plox requires a quadratic root to be found, M T ppl °* (like Mct) does not require iteration to 
compute, and is found to perform better then Mct in circumstances in which the information from 
the missing transverse momentum (which the former retains and the latter discards) is both reliable 
and useful. 



I. INTRODUCTION 



The stransverse mass, -Mt2, [I] is a kinematic variable which is defined to be the maximal lower-bound on the mass 
of a parent particle decaying to an invisible daughter of a known (or hypothesised) mass and one or more visible 
daughter particles in a hadron collider event, provided that there are two of these parents in the event. Properties 
and generalisations of Mt2 have been investigated extensively [2 25 . Experimentally, Mti has been used in searches 
for R-parity conserving supersymmetry signals [26 -30J, Higgs searches (for example [31, 32 ) and in measurements of 
the top quark mass |33) . 

Although its properties as a mass-scale identifier and QCD-background discriminator are well established, it is not 
possible in general to write down a closed-form analytic expression for Mt2 or many of the variables to which it is 
related. For some discussion on this point see 0[17l|34]. To date therefore, Mti and similar variables are generically 
evaluated using iterative numerical techniques, and code is available from various sources to do this [3SJ EH]- It has, 
however, been shown that in a wide class of special cases [37] analytic formulae for Mt2 can be found. In one of these, 
specifically the so-called "fully massless case" in which the daughter particle masses are all assumed to be zero, Mt2 
may be written as a simple function of one of the two or four real roots of a quartic equation: 

As 4 + Bs 3 + Cs 2 +Ds + E = (1) 

whose coefficients are simple functions of the observed momenta. 1 On its own, this result is perhaps little more than 
a curiosity. Quartic roots are not terribly pleasant things to work with - so much so that they are frequently treated 
numerically even though analytic formulae for them exist. As a result, very few implementations have made use of 
analytic solutions of equation given fears of instability, and it is even less clear that any useful insights are likely 
to emerge from pursuing this line further. The idea that motivated this investigation was a desire to challenge that 
assumption — a desire to see if anything useful was trapped inside equation (fTl) after all. 



'Electronic address: lallyc@cranbrook.kent.sch.uk 
^Electronic address: lester@hep.phy.cam.ac.uk 

1 Explicit expressions for the hve coefficients appear later in the paper. 
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FIG. 1: The diagram shows an initial state (the left most blob) which produces a pair of parents (A and A') and some other 
source (or sources) of Upstream Transverse Momentum (g). Each visible particle decays to one or more visible particles (B) 
and an invisible particle (C). 

The starting point, therefore, requires the construction and validation of a stable method of calculating the roots 
of equation ([T]) using non-iterative analytic methods (described in detail in the appendix). The important part of 
that task is not the root finding itself, 2 but rather the understanding of the characteristics of the input space that 
determine the characteristics of the roots - and in particular determine which of those roots is the "physical" one that 
is ultimately used to calculate Mt2- 

This paper documents the insights that were gained regarding the geometric properties of M^2 found in the special 
case considered, and indicates how they could used to define an approximation to Mt2, computable by a fixed-cost 
non-iterative method, possibly as an alternative to Mqt appropriate when the missing transverse momentum is 
reliable input, and where fast and reliable calculation times are required, for example in trigger routines. 

II. REVISITING THE MASSLESS CASE 

Following the notation of [37] we revisit the derivation of the expression for Mt2 in the massless case that uses ([!]) , 
with a view to highlighting the parts that are related to how the correct choices are made among those roots. 

We begin with the hypothesis, illustrated in figure [l] that two pair-produced parent particles (A, A') each decay 
into a (group of) visible (B, B') and invisible (C, C) daughter particle(s). Each decay branch is referred to as one 
"side" of the event. The momenta of the visible decay products of sides 1 and 2 (assumed to have been reconstructed 
into two jets) are referred to as a M and 6 M respectively, and the hypothesised momenta of the invisible decay products 
as p M and respectively. 

We restrict ourselves to the transverse plane and so we view a M as consisting of three components: a M = (e-TiPxiPy), 
or = (ct, a). Furthermore, as we assume a fully massless situation, then a 2 (= e\ — a 2 ) = b 2 = p 2 = q 2 = 0. 

Of course the invisible momenta p^ and q^ will not be known; instead for each event the missing transverse 
momentum jzi will be obtained by adding together the visible transverse momenta, a and b, and the "upstream 
transverse momenta" (UTM), denoted by g. Thus we have the following relationships between the transverse momenta: 

a + b + j* + g = (2) 

p + q = (3) 

Equation |3j is referred to as the "momentum splitting condition" . Without loss of generality, we assume that lies 
outside of the sector of the transverse plane spanned by positive multiples of a and b. 3 We also introduce the (2-D) 



2 Nonetheless, ensuring stability in any root-finding code is non-trivial, and the appendix also details why the method employed here is 
regarded as being sufficiently stable. 

3 We remind the reader that in cases in which the jzt vector happens to lie "between" a and b then Mj-2 i n the "fully massless" case must be 
identically zero (the so-called "trivial zero" of Mt2) as we can satisfy the momentum splitting condition by taking p and q collincar with 
a and b respectively. This trivial case, in which Mt2 is identically zero, is implicitely excluded from any of the subsequent discussion 
concerning methods of evaulating Mt2, though it is always present in any actual evaluation or implementation. (See right-hand plot of 
figure [6| later for an example of this "trivial zero" scenario.) 



3 



unit vectors of the transverse momenta (actually these are the transverse velocities of the particles as in the massless 
case e T = |p|): a, b, p, q, p. 

Mt2 itself is defined as the minimum (over all hypothesised momenta p and q consistent with the momentum 
splitting condition) of the larger of the two transverse masses obtained from each "side" of the event. Symbolically: 

M T2 = min {max [M T (a M , p) , M T (b^ , q)] } (4) 
p+q=P» 

For more details see pQ. For now all we need to note is that this is a min-max process, and in calculating Mt2 our 
interest will be drawn primarily to the specific vectors p and q which minimise the max in the RHS of eq Q. We 
will refer to these vectors as those which "realise the extremum of the Mt2 min-max process" , or sometimes simply 
as those vectors which "lead to the Mt2 solution" etc. 

A key property of the hypothesised vectors q and p that realise the extremum of the Mt2 min-max process, is that 
the angle between them is the same as the angle between a and b [37 . As a consequence, p and q can be parameterised 
in terms of a single unknown parameter in the following way: 

p = cos Oh + sin Oh' 

q = cos Ok + sin Ok' (5) 

where the two vectors a' and b' are defined to be the two vectors a and b rotated by +90 degrees in the transverse 
plane. In effect, represents the angle by which a and b must be rotated by in order to arrive at the directions of q 
and p which lead to the Mt2 solution. 

As we have assumed that p lies outside of the sector of the transverse plane spanned by positive multiples of a and 
b, we are only interested in the so-called "balanced" momentum-hypothesis configurations where both sides of the 
event have the same transverse mass. In other words, p and q must also satisfy the following: 

2(|a||p|-a.p) = 2(|b||q|-b.q), (6) 

and when this contraint is met, each side of equation will also be equal to the desired final value of M|, 2 . Putting 
the "opening-angle" condition ^ together with the balanced hypothesis" condition (|6| leads to the first of two 
useful results: 

|a| |p| (l - cos 6»(a.b) - sin 0(a.b')) = |b| |q| (l - cos H^) + sin 0(a.h')) (7) 

while taking the opening-angle condition ([5| together with the momentum splitting condition ([3| and a dot-product 
with p' (p 1 rotated by +90 degrees in the transverse plane) leads to the second: 

+ |p| (cos0(h.p) +sin0(b.j*)) = - |q| (cos 0(a.j*') + sin0(a$)) . (8) 

The quotient of these two constraints Q and {[sj) eliminates |p| and |q| allowing to be found. Indeed, it is 
precisely that quotient re-written as a quartic in terms of s = sin 0, that was referred to in equation ([I]) wherein, for 
completeness, we note that 

A = Kj + (K ss - K cc f (9) 

B = 2(K CS K C + K S (K SS - K cc )) (10) 

C = K s 2 - Kj + K C 2 + 2(K SS ~ K CC ){K X +K CC ) (11) 

D = 2(-K cs K c + K s (K 1 +K cc )) (12) 

E = (Ki - K c + K CC )(K 1 + K c + K cc ) (13) 



using the definitions 



K ss = -(S.p)(a.h') (14) 

K cc = -(a.p)(a.b) (15) 

K s - {a.p) (16) 

K c = (a.p ) (17) 

K cs = -(a.p)(k.h) - (S.p)(k.h') (18) 

K x = (19) 

a = a + b and (20) 

5 = a-b. (21) 
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(a)Typical event with two real roots 



(b)Typical event with four real roots 



FIG. 2: Unit-vector representation of momenta a, b and p, and candidate p and q values obtained from the roots of the 
quartic 



As before, momenta unit-vectors carry a "hat" whereas a "prime" (as in b') signifies a rotation of the unit-vector 
through +90 degrees in the transverse plane. 

The key to understanding which of the real roots (i.e. always either two or four real values for sm9) leads to 
the correct physical value for Mt2, is this geometric picture of what we are doing to calculate Mt2- rotating the 
hypothetical invisible p and q momenta away from the visible momenta a and b until both the "balanced hypothesis" 
constraint Q and the momentum splitting constraint ^ are satisfied. What in effect this means is that the correct 
root must correspond to an angle that rotates p and q until one and only one of these vectors has been rotated "past" 
p' thus ensuring constraint (|8| can be satisfied. 4 Most of the time this rotation condition is sufficient to identify 
the correct physical root as most of the time we only have two real roots (the other two complex roots having been 
discarded as obviously unphysical) . Figure 2(a)| taken from a toy Monte Carlo event generator used to test the quartic 
root finder code (described later) , shows a typical example of this situation where it is clear which of the two possible 
rotated pairs of vectors (i.e. the pair numbered "1" or "2") would satisfy the momentum splitting condition |3| 5 . 

There are occasions where we have four real roots and more than one satisfies the rotation condition above. These 
all potentially correspond to physical values of M T2 . But as M T2 by definition requires a lower-bound of the original 
parent masses, the final condition we require is that the correct root is the one whose corresponding angle 9 ensures 
p and q are as parallel to a and b (respectively) as possible, as this will identify the lowest value for Mti for the 
event in question (i.e. this by definition minimises both sides of equations (|6| and ([7|). Figure 2(b) shows a typical 
example of this situation, and here, although vector pairs "1" , "2" and "3" correspond to possible physical rotations 
of p and q 6 , it is vector pair "3", i.e. vectors P3 and q3, which are most "parallel" to a and b and so it is this pair 
which corresponds to the lowest Mt2 value. Note that in very many cases, this is the same as saying it is the root 
whose rotation angle is as close to 180 degrees as possible - more on this point later. 



4 Note the rotation is defined to be always in the "positive" , that is anticlockwise, sense. 

5 It is worth noting that the other real root (e.g. in this case the vector pair numbered "2" in figure [2~( a) | corresponds to values for |p| and 
|q| which have negative values. These still satisfy both the balanced hypothesis constraint ([TJ and the momentum splitting constraint 

f .|, but are clearly unphysical. 
ut not vector pair 4 as }S does not lie "in-between" them and so they represent the unphysical solution for 6. 
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FIG. 3: Low UTM scenario input momenta and output rotation angles 



III. AN IMPLEMENTATION 



Using a method due originally to Leonard Euler - but expanded upon recently by Nickalls |38] whose work was 
especially useful in ensuring a stable root-finding method (detailed in the appendix) - and the right root identification 
techniques outlined in the previous section, a C++ program was written to calculate Mt2 hi the fully massless case 
(the code can be found at [35]). 



Testing and Accuracy 



In order to test this implementation, a simple "toy" Monte Carlo was used to generate signal-type events for the 
program to analyse. The toy Monte Carlo reproduces a simple event as depicted in figure [T] it begins with a normally 
distributed centre-of-mass "rest mass" which converts isotropically into two pair-produced parent particles (each 
possessing half the rest mass of the original "rest mass" ) , which themselves decay into daughter visible and invisible 
particles (the latter's energy and momenta being represented by the missing energy and momenta of the event). For 
the purposes of this analysis, the mean of the normally-distributed pair-produced parent masses was fixed at 750 GeV, 
and the visible and invisible decay particles were assumed to be massless. This leads to events with a wide range 
of momenta for the visible and invisible decay particles. UTM in an event is introduced by including a normally 
distributed initial centre-of-mass energy (which can boost the momenta of the visible and invisible decay particles); 
adjusting this initial centre-of-mass energy distribution up or down produces events with more or less UTM. 

The initial "low-UTM" run uses a mean centre-of-mass energy of 200 GeV which produces UTM magnitudes that 
are on average 50% of the size of the missing transverse momentum in an event. Figure 3 (a) | shows the distributions 
of visible daughter transverse momenta, missing transverse momenta and event UTM, for a typical 10,000 event 
"low-UTM" run. 

The algorithm was tested against many hundreds of thousands of signal-type Monte Carlo generated events of the 
form described above, and has shown good agreement with the best numerical techniques available: the fraction of 
evaluations in which the algorithm was found to generate Mj-2 values differing by more than O.OlGeV (or by order 
0.01%) from the known correct values for the events generated (see later description) was found to be O(10 -5 - 10~ 4 ). 



Observations 



It is interesting to note that the root values (remember these are actually equal to sin 9 where 9 is the angle required 
to rotate vectors p and q away from a and b until equations ^ and ^ are satisfied) are very skewed towards a value 
for 9 of 180 degrees. Figure 3(b) shows the "low UTM" distribution for cos 9 (cosine of the rotation angle is used for 
clarity) demonstrating that an angle of approximately 180 degrees (i.e. cos 6* = —1) is highly favoured. 

To understand why ~ 180 degrees is the dominant solution to the Mti quartic expression in the fully massless case 
we need to look again at the geometric picture of what is happening with respect to momentum vectors a, b and p 
(i.e. p + q) depending on the exact size and direction of the UTM of the event. 
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FIG. 4: Two sample events showing the transverse momenta unit-vectors where UTM is negligible. The left-hand plot shows 
an event where p is relatively large compared with the visible momenta a and b (thus a and b are "pushed close together"). The 
right-hand plot shows an event where conversely is relatively small versus a and b (thus a and b are more "back-to-back"). 
The unit-vectors for the invisible momenta p and q, determined by evaluating the quartic expression for the events, are also 
shown. Both of these events have a value of 180 degrees for the rotation angle 9. 

A. An Important Result when the UTM = 

An important configuration occurs in this massless case when the UTM is zero (or negligible) . Equation ^ becomes 

a + b + p 1 = (22) 

hence 

a + b = -p-q (23) 

In order to continue to satisfy the requirement that the angle between q and p for the Mt2 solution is fixed to be the 
same angular separation between a and b (and given the parameterisation of (J5|) we have in this situation that 



a = 



= |q| and |b| = |p| 



(24) 



and that 8 is exactly 180 degrees i.e. in the fully massless case, when UTM is zero, in order to conserve momentum, 
the parent particles decay such that the visible daughter of one "side" of the event has an equal but opposite momentum 
to the invisible daughter from the other "side" of the event. This means that the value of Mt2 in the massless case 
when UTM is negligible, given by either side of expression 0, simplifies to 



2|a| |b| (l - cos 0(a.b) - sin6>(a.b') 



(25) 



hence 



2|a||b| l + (a.b) 



(26) 



M| 2 



2|a||b|(l 



COS ( 



(27) 
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FIG. 5: Two situations where large UTM has (left plot) little effect, and (right plot) a bigger effect, on the rotation angle 
leading to the Mt2 solution. 



where <fi is the angle between vectors a and b. The expression on the right-hand side of this equation is actually also 
that for the so-called contransverse mass Mq T as defined in |12) 7 . 

Figure [4] shows some sample Monte Carlo events where the UTM is negligible. It is clear visually why conservation 
of momentum leads to the relationships of (23) and (24 1 between a, b, p and q. 



B. When the UTM > 



The interesting result presented here showing that a majority of events tend to have 9 close to 180 degrees is due 
to the fact that even when the UTM is non-zero, it will not necessarily change the rotation angle significantly. This 
is because it is not just the magnitude, but also the specific direction of the UTM w.r.t. the directions of the visible 
and invisible momenta that determine if the UTM will have any affect on 9. Specifically, the UTM would need to be 
both sufficiently non-parallel to one (or more) of the daughter particle momenta, but also of the order of or larger 
than that specific momentum vector's magnitude. Visually this can be thought of as the UTM "pushing" a smaller 
momentum vector to one side. 

For example, in figure[5j the left-hand plot is very similar to the left-hand plot in figur e [4| (i.e. large p and relatively 
small visible momenta), except that here the UTM is significant (160 GeV versus 190 GeV for p, 48 GeV for a and 
58 GeV for b). But because the UTM is roughly anti-parallel with p, it has virtually no effect on the direction of 
the missing transverse momentum and a minor effect on the angular displacements of a and b (actually the visible 
momenta will have been "pushed away" from each other somewhat compared with if there had been negligible UTM) , 
and thus the appropriate rotation for p and q remains at close to 180 degrees. For this type of event to have a rotation 
angle significantly different from 180 degrees, the UTM vector would have needed to be at a substantial angle to the 
other momenta vectors i.e. in terms of momentum phase-space, only a relatively moderate sub-set of vector directions 
of UTM would be able to achieve this. The right-hand plot of figure [5] shows an event of this type where the UTM 
has "pushed" vector b (and also in this case the missing transverse momentum) away from it, thus decreasing the 
rotation angle for p and q to around 138 degrees. 

One final comment to make on how UTM may or may not affect the rotation angle is to look at situations similar 



7 This equality between M^2 an d A^CT i n the fully massless case when the UTM is zero has already been noted by Serna 1141 using a 
more general derivation. 
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to the right-hand plot in figure H (i.e. where p> is relatively small versus a and b and thus a and b are roughly 
"back-to-back" ) . The left-hand plot of figure pi shows again that a significant UTM is not on its own sufficient to 
move the rotation angle away from 180 degrees. Here a UTM of magnitude 71 GeV is roughly (anti-) parallel to the 
back-to-back a and b vectors (101 GeV and 36 GeV respectively), has therefore virtually no effect on their directions, 
and thus has no affect on the resulting rotation angle which continues to be approximately 180 degrees. 

However in this type of back-to-back event, even if the UTM has the right direction to "push" a, b or p towards 
each other, if this were to cause to now lie between a and b, then we are back to the trivial zero situation where the 
balanced value for Mt2 is identically zero on both sides of the event as mentioned in section [Tl] (the right-hand plot of 
figure [6] shows this type of event). So in this low ]/>, back-to-back situation, there is an even more limited phase-space 
for the UTM vector to move within in order to produce a rotation angle significantly different from 180 degrees. 8 



IV. APPROXIMATING M T2 USING A FIRST ORDER CORRECTION TERM TO 6 



Having understood that in the "fully massless" case there is a strong bias towards the rotation angle 9 being close 
to 180 degrees, it is sensible to consider if the angle which leads to the correct balanced value for Mti could be 
approximated by assuming 9 is equal to 180 degrees plus a first order correction term; this approximation might 
possibly circumvent the need to solve the quartic equation at all. 

To investigate this we need to look again at the constraint equations Q and ^ whose quotient eliminated |p| and 
|q| and which lead to the quartic expression. Looking at ^ when 9 is 180 degrees (remembering that this signifies 
an event where UTM is negligible) it is clear that it is the vanishing of the sin 9 term that allows the expression to 
be simplified (and led us to the derivation of the Mqt equivalent expression of (27l). When 9 is not 180 degrees, it 
will therefore be this, now non-zero, sin (9 term which, as it has the opposite sign on either side of ([7]), will have the 
greatest effect on adjusting ^ (and also (j8J)) to allow the correct balanced values of Mt2 to be calculated and the 
equality to continue to be satisfied. 

As an approximation, we will therefore assume that the rotation angle is represented by (9+59), where 59 is a first 
order correction term for the angle. In what follows we will refer to this estimated rotation angle as 9 e (and to the true 
rotation angle, i.e. as determined from the quartic roots, as 9 t ) and we assume as a starting point that 9 e equals 180 



Note therefore that in situations of higher UTM, we would expect there to be more trivial zero type events (i.e. Mt2 values of zero). 
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degrees. As 89 is by definition a small angular adjustment, then to a good approximation the cos(6> + 89) terms can 
be replaced by -1 and the sin(# + 89) terms can be replaced with -89. The two constraint equations would therefore 
simplify to: 

|a| |p| (l + (a.b) + <J0(a.b')) = |b| |q| (l + (a.b) - <J0(a.b')) ( 28 ) 



+ |p| {-(b.p ) - S9(b.f)) = - |q| (-(a.p ) - <W(a.j*)J (29) 

Taking the quotient of both sides (again to eliminate any dependence on |p| and |q|) and rearranging we obtain a 
quadratic function in 89: 

B(\b\D - \a\F)89 2 + (\a\(BE - AF) - \b\(AD + BC))89 + A(\a\E + \h\C) = (30) 

in which 

A = (l + (a.b)) (31) 

B = (a.b') (32) 

C = -(£.£') (33) 

D = (b.f) (34) 

E = -(a.p) (35) 

F = (36) 

This quadratic equation can be solved using the usual quadratic formula 9 and the value for 89 thus gives a correction 
to our assumed 180 degree value for 9 e . This corrected value can then be used to calculate approximate balanced 
values for Mt2- We will refer to these approximated values going forward as M^ proK values. We would expect in 
the cases where 9t is close to 180 degrees, that our approximated value for 9 e , and thus M^ prox , would be very close 
to the "true" values calculated using the quartic algorithm. In the rarer events where the UTM produces a rotation 
angle significantly different from 180 degrees (say more than 45 degrees different), we would not expect the correction 
term to produce a particularly accurate estimation of 9 t and thus Mt2- However in these events, we can again use 
the geometry of the massless case to estimate the correct rotation angle: 



As discussed in section [TTTJ it is the symmetric distribution of the visible and invisible momenta in low UTM events 
which leads to rotation angles close to 180 degrees. Conversely because events with 9 t significantly different from 180 
degrees is caused by a UTM vector which "pushes" the visible and invisible momenta closer together, these events 
are by definition asymmetric and tend to have a and b vectors "close together" (the right hand plot of figure [5] 
shows this typical situation) . We can make a simple estimation of the correct rotation angle by choosing one of three 
possible scenarios (i) a rotation angle that is just more than the minimum possible angle (i.e. that which ensures that 

is just inside p and q), (ii) a rotation angle which is just less than the maximum possible angle (again ensuring 
p 1 remains inside p and q), or (iii) a rotation angle that sits halfway between p and q. To determine which is 
most appropriate we can, for example, look at the relative magnitudes of |a| and |b|, as for the balanced situation 
(i.e. see equation ({61)) we would expect a dominant a vector (say where |a| is more than twice the size of |b|) to be 
associated with a dominant q (i.e. from the parametrisation of q in (JHJ) and a dominant b to result in a dominant 
p (from equation ([5])); the right hand plot of figure [5] shows an example of this latter scenario. By looking at the 
relative magnitudes of |a| and |b| for an event, we can therefore choose the minimum or maximum angle that ensures 
the appropriate vector dominates 0, and thus estimate Mt2 for that event. If neither visible momentum vector is 
dominant, we choose scenario (iii) which estimates equal |p| and |q| values. 

By adding this estimation method to the M^2 Pr ° x approximation, we now have an algorithm which we hope can 
provide a fast calculator for estimating Mt2 irrespective of the rotation angle in question. In order to check the 
algorithm's accuracy, in addition to testing it in the "low UTM" scenario discussed previously, the same toy Monte 



The sign of the 88 term is assumed to signify which of the two roots is the correct physical value (i.e. whether to choose the + or — in 
the quadratic formula). 
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"High UTM" Transverse Momenta Distributions 
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(a)Momenta distributions in the "Mid UTM" scenario 



(b)Momcnta distributions in the "High UTM" scenario 



FIG. 7: Momenta distributions of higher UTM scenarios 
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FIG. 8: cos 8t distributions of higher UTM scenarios 



Carlo generator was used to produce events in two other scenarios: the "mid UTM" run where the UTM in an event is 
on average equal to the typical p magnitude, and the "high UTM" run, where the UTM is on average twice the size of 
the typical ^ magnitude. The distributions for visible daughter transverse momenta, missing transverse momenta and 
UTM for these higher UTM scenarios are shown in figure [7J The resulting values of cos 9t obtained for these two runs 
using the quartic calculator are shown in figure [8] The longer tails should be compared with the previous figure |3(b)| 
demonstrating that scenarios of higher UTM are more likely to produce rotation angles significantly different from 
180 degrees. However if should also be noted that even in these higher UTM regimes, the rotation angle distribution 



is still dominated by events with 9 t close to 180 degrees as expected and as discussed in section III 



The approximate method was used to calculate M^ pTOK for 10,000 "mid UTM" and "high UTM" events as well as 
the original 10,000 lower UTM events, and these were compared with the "true" values as obtained with the quartic 
calculator. Firstly Figure [9] shows the resulting distribution of the percentage difference between the "low UTM" 
^appiox va i ue anc i "t rue " Mt2 value for each event (percentage difference calculated by subtracting the "true" 
Mt2 value from the M^ prox value and dividing by the "true" M^i value - number of events per bin plotted on a log 
scale to see the distribution tail more clearly) . The accuracy of the Af^2 Prox algorithm is very encouraging with over 
95% of events within 2% of the "true" Mt2 value. 
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FIG. 9: Percentage difference between M^ vo * and Mti in the "low UTM" scenario 
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FIG. 10: Percentage difference between Myj 11 ™ and "true" Mt2 in the higher UTM scenarios 



Figure 10 shows the same comparisons for the two higher UTM runs 10 . As would be expected from the higher 
UTM in these events (leading to more events with rotation angles significantly different from 180 degrees), the M^2 Prox 
values have larger tails than the "low UTM" situation but still demonstrate very good accuracy: for the "mid UTM" 
run about 90% of values lie within 2% of the "true" Mt2 value, and for the "high UTM" run over 70% of them are 
within 2% of the "true" value (and 90% within 10% of the "true" value). 

Although the MtjSf* * calculator demonstrates good accuracy for estimating the Mti value for an event, the real 
benefit it brings is one of speed. 



V. COMPARISON OF THE SPEED OF THE "TRUE" M T2 AND A/^ 2 PPROX CALCULATORS WITH THE 

FASTEST NUMERICAL CALCULATOR 

For the "true" My 2 calculator, where Mti is determined from the physical root of the underlying quartic equation, a 
typical calculation (where transverse momentum values a, b and ^ are provided to the code) takes a little less than 3/is 
per event to run on a standard PC (ignoring any read/write times to and from datafiles). This is approximately 4 times 



10 In the "high UTM" run, it is not surprising that approximately 16% (1,626) were of the "trivial zero" type with Mj>2 = (compared 
with less than 10% (837) in the "mid UTM" run). The histograms here therefore exclude these "trivial zero" events and show 8,374 
events (high UTM) and 9,163 events (mid UTM). 
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FIG. 11: Comparison of M™ prox and Mct as mass-scale identifiers 



faster than the current fastest numerical Mt2 calculating code provided by [34] in their ellipse bisecting procedure 
(used in the massless case). One limitation on speed for the "true" Mti algorithm seems to be the trigonometric 
nature of the quartic solving procedure requiring a number of library calls within the code. The authors are inclined 
to believe that more experienced programmers could modify the current code and improve this initial performance. 

For the Al!^ plox algorithm, however, this analytic method can calculate a good estimate of Mt2 m l ess than 1/xs 
per event on a standard PC; that is more than 10 times faster than the current fastest numerical calculator. 



VI. USE OF M T 2 IN THE FULLY MASSLESS CASE AS A TRIGGER VARIABLE 



As discussed in |37j , one of the motivations for finding useful "special cases" in which the solution for Mt2 is fully 
analytic and can therefore be evaluated quickly and accurately, is its possible use as a "trigger" variable at the LHC. 
The properties of Mt2 as a variable which is sensitive to the mass-scale of the pair-produced parents and its ability 
to discriminate larger mass-scale events from low multiplicity, relatively low-mass-scale, QCD events, would indicate 
that Mt2 should perform well as a trigger to help the LHC experiments cope with the ever-increasing instantaneous 
luminosity conditions. 

As has been mentioned, the new analytic calculator presented here is the fastest yet for calculating Mt2 in the 
fully massless case (and improvements in speed may be possible still), and with regards to the M^ lox approximation 



discussed in section |rV] this has been shown to be a good approximation to the true Mt2 value and, as a much simpler 
calculation, may be fast enough to compete with other possible trigger variables available. 

As an example of its sensitivity, we have compared its performance as a mass-scale indicator against that of Mct 
(the contransverse mass, which in section III we noted is equal to Mt2 in the fully massless case when the UTM is 
exactly zero). For the "mid" and "high UTM" runs, we compare the value of Mct and M^2 P1 ° X against the rest 
mass of the pair-produced parent particle giving rise to the decay being analysed (remember this was chosen to be a 
normally distributed set of mass values with a mean of 750 GeV). Figure 
M CT and M^ p 2 plox and the parent mass for both the "mid UTM" and "higFUTM 
its construction as a maximal lower-bound, Mt% would never give a higher mass than that of the parent and it is 
pleasing to see that M™ pr ° x also demonstrates this quality, all events being less than or equal to zero mass difference. 
However Mct has a significant tail demonstrating that, based as it is on the visible momenta of an event only, that 
when UTM is non-zero and thus the magnitudes of the visible momenta are larger than they otherwise would be to 
balance this UTM, Mct tends to overestimate the mass-scale of an event. 



shows the mass difference between both 
runs. As would be expected from 



VII. CONCLUSIONS 



We have constructed a fast and stable analytic calculator for Mt2 in the "fully massless" case (the scenario in 
which Mt2 has been most frequently used at the LHC) based on the solving of a quartic equation. Unfortunately 
due to the trigonometric nature of the quartic solution and the necessary library calls required by the code, the speed 
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of this calculator is only around 4 times faster than the fastest available numerical code due to [53] (though we will 
continue to work on the algorithm and believe that further improvements in speed are achievable). However due to 
the kinematics characteristic of the "fully massless" case we have also been able to approximate Mt% using a much 
faster and more straightforward algorithm that requires little more than the solving of a quadratic equation. This 
approximation, M^™, has been shown to be an accurate estimator of Mt2, and Mj^ values can be calculated 
more than 10 times faster than the current fastest numerical calculator. 

It was the authors' hope in beginning this analysis of M^ 2 in the "fully massless" case that it might lead to a 
Mt2 calculator that could be fast enough to use as a "trigger" variable at the LHC and next linear collider. The 
^appiox gjgQj-^jjjQ ma y De a possible step in that direction and it's usefullness as a mass-scale identifier has been 
demonstrated. 
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IX. APPENDIX : QUARTIC ROOT FINDING 

Methods for solving quartics analytically 11 have been known for centuries and Ferrari is credited with achieving it 
first in 1540 [39]. The methods (all essentially based on Ferrari's original [40]) involve manipulation of the quartic 
coefficients in order to end up with an easier to solve cubic expression called the "resolvent" cubic. As with all 
polynomials, the roots of cubic and quartic equations can be real or complex. In fact the three roots of the resolvent 
cubic are either all real (leading either to four real quartic roots or to two pairs of complex conjugate quartic roots), 
or have one real and two complex conjugate roots (leading to two real and one pair of complex conjugate quartic 
roots). 

Unfortunately, due in part to the number of algebraic steps required, algorithms for solving quartic equations tend 
to suffer from instability [41] . These instabilities can arise due to, for example, numerical round-off errors in the 
calculations or the handling of imaginary numbers. Two situations relevant to this analysis are: 

• where the resolvent cubic has three real positive roots (and the quartic should therefore have four real roots) 
but rounding errors cause one or more of the smallest of these cubic roots to appear to be negative - this error 
will lead the algorithm to obtain incorrect complex roots for the quartic. 

• where the resolvent cubic has one real and two complex conjugate roots but rounding errors cause the imaginary 
values of the complex conjugate pair to be unequal - this will lead to four complex roots for the quartic, when 
in fact two of the roots should be real. 

As we know a priori that the quartic expression for Mt2 must have at least two real roots (i.e. one of the two real 
roots will lead to the physical value for Mti), we require a method that makes explicit (a) whether or not a given 
quartic has any complex roots - this enables us to identify immediately rounding errors such as in (i) above 12 , and 
(b) the process by which the three resolvent cubic roots are combined to produce the four quartic roots - this allows 
us both to correct any rounding errors to the imaginary parts of the resolvent cubic's complex roots, and to minimise 
the need to manipulate the imaginary parts of any complex roots, as we as we know these must cancel when forming 
the two real quartic roots (i.e. we are only interested in any real roots). 

A method developed by Leonard Euler in the eighteenth century [4121 1415] , recently expanded upon by Nickalls [55J , 
allows the following: 

• A simple way to calculate the discriminant of the quartic equation thus showing immediately if the roots of the 
quartic are all real or if there are two real/two complex roots (and thus also determining if the resolvent cubic's 
roots are all real or one is real and the others form a complex conjugate pair). 



Quartics being the highest degree polynomial which can be solved analytically. 

i.e. if the quartic has all real roots then we know the resolvent cubic must have three real positive roots and we can proceed accordingly. 
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• A relationship between the resolvent cubic roots and those of the quartic is given explicitly thus allowing a 
stable algorithm to be constructed that derives the latter from the former without concern for rounding errors 
in, or any significant manipulation of, imaginary numbers. 

Furthermore, Nickalls |44j gives two complementary methods for calculating the roots of the resolvent cubic de- 
pending on whether or not there are one or three real roots. The benefit this brings to the algorithm is that in the 
situation where we know we have only real roots, we can use a more stable trigonometric calculation to obtain them. 

The rest of this appendix summarises the approaches detailed in [38] and [44] and which underlie the resulting code 
for calculating the roots to the M^ 2 quartic expression. 



A. Resolvent cubic 



Leonard Euler was the first person to fully appreciate to role of the so-called "resolvent" cubic in calculating the 
roots of quartic equations [32 03]. Nickalls has written at length on the geometry of, and solutions to, various 
polynomials, and the details in this appendix are from two particular papers [38j S3]. The notation follows that in 
[58] , but sometimes also shows the notation in [33] for ease of reference 13 . 

For any given quartic of the form 

ax 4 + bx 3 + cx 2 + dx + e = (37) 

manipulation of the coefficients of the quartic allows the "reduced" form of the quartic to be derived (the reduced 
form has no x 3 term). This is easily done by putting x — y + A/v/, where Ajv f is equal to (— 6/4a), and substituting 
for x in the original quartic expression. Euler's quartic solution showed that each of the roots of the reduced quartic 
is the sum of three square roots, ±^/r\ ± ± y^rjj, where t\ etc. can be derived directly from the three roots of an 
underlying "resolvent" cubic (which we will call 01,02,03). A convenient form of the resolvent cubic is the reduced 
form (no x 2 term) as follows: 

T{z) = z 3 - ZIz + J = (38) 

where the quartic invariants, / and J, are given by 

I = I2ae-3bd + c 2 (39) 
J = 72ace + 9bcd - 27ad 2 - 27eb 2 - 2c 3 (40) 

Furthermore, the invariants I and J can be combined to form the discriminant of the quartic equation (4/ 3 — J 2 ) 
and thus determine whether or not the resolvent cubic has three real roots (when (4/ 3 — J 2 ) > 0) or one real root and 
one pair of complex conjugate roots (when (4/ 3 — J 2 ) < 0). In the quartic calculator, we use two different methods 
for finding n etc. depending on whether not the resolvent cubic has only real or some complex roots. 



B. Real Roots 



If the resolvent cubic has three real roots (which we know must lead to a solution for our quartic equation of four 
real roots as one of those roots must lead to the real, physical value for Mx2), then we can use a trigonometric method 
for calculating the roots as follows (as per [33]) : we define an angle, 8, such that 

cos 30 = — — = ; (41 

h VI7 3 

Then the three roots to the reduced resolvent cubic are given by: 

ci = 2(5cos<9 (42) 

c 2 = 25 cos (43) 

in 
3~ 



c 3 = 26 cos (44) 



The rationale for, and derivations of, expressions which appear in this appendix can be found in the relevant references. This appendix 
will simply set out the procedure for finding the quartic roots. 
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where S = \fl. The Euler values r±, T2, ^3, are then equal to: 

r <=' 2 + k < 45 > 

+ < 46 > 

where e 2 = — ^-5— — . Note that as we know the quartic expression that leads to Mt2 must have at least two real 
roots, we know that n, ri and r% must all be positive numbers (i.e. if any were negative, this would lead to four 
complex quartic roots). 



C. Real and Complex Roots 



If the resolvent cubic has one real root and one pair of complex conjugate roots then we must use a method which 
can handle complex numbers. The real root is given by: 



Cl 



and the pair of complex conjugate roots are given by: 



C2 



C3 



2 

ci 
" 2 



2 

.V3 



a/ci 2 - 4/ 



4/ 



4/3 



(48) 

(49) 
(50) 



The three Euler values r%, ra, r3, are then also given by (45), (461 and (|47| (but now we must remember that 
and r"3 are complex numbers). 



D. The Quartic Roots 

As has been mentioned previously, the four quartic roots are obtained from four of the eight possible permutations 
of iy^i y/r^ ± A/F3. Euler showed that for a given quartic equation, as the product of y^y^y^ would either be 

a positive or negative value, that the correct sign would be the same as the sign of (-^-) where Y Nf is given by 
the following expression: 

b 3 +8a 2 d~4abc 

Ysf> = ^ (51) 

Thus once we calculate the sign of ^ Y " f ' 1 from the quartic coefficients, we know which four permutations arc the 

ones that lead to the correct roots of the original quartic equation ( |37| (which we will call X\, X2, X3, X4). So we 
have the following two possible results: 

lf(z^)<0 

(52) 
(53) 
(54) 
(55) 



And when > 



Xi 


= X N f + 


s/n + 


\/^2 - 


A3 


X 2 


= x Nf + 


\fn - 


s/T2 + 




x. 


= X N f - 


x/n + 


s/^2 + 




x 4 


= X Nf - 


vAT- 


\pT2~ - 


V^3 


X! 


= X Nf - 






A3 


x 2 


= X Nf - 


s/n - 


vA^ + 




x 3 


= Xpff + 


s/rl - 






X A 


= Xjff + 


v/n~ + 


\/r2 + 


V^3 



(56) 
(57) 
(58) 
(59) 
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Finally we note that in the case where we know that the above must give two real and one complex conju gate pair 
(i.e. when (4/ 3 - J 2 ) < 0), then we know that the two real roots we want must be derived either from rt54| and 
( 55 1 [or from ( 58 ) and ( 59 ) as appropriate] because only in these expressions will the imaginary parts of ^Jr^ and 
-y/?"3 cancel. Thus at this stage not only can we ignore (52) and (53) [or (56) and 1571)], but we can also ignore the 
imaginary parts of (54) and (55) [or (58) and (59)] because these must cancel. This is not just expedient, it also 



minimises any stability issues in the a 
numbers (as discussed previously). 

Of course if we know we instead have four real roots (i.e. when (4/ 3 



gorithm that might arise from errors in accurately calculating these complex 



J 2 ) > 0), then we must evaluate all four 



roots as any one of them might be the one that leads to the minimum value for Mt2 ■ 
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